Competition between spin glass order and strong coupling superconductivity in a 

single - species fermion model 
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The phase diagram of a single species fermion model allowing for local pairing superconductivity 
(SC) and spin glass order (SG) is derived as a function of chemical potential /x and ratio r = 
v/J between attractive coupling v and frustrated magnetic interaction J. For ratios larger than a 
characteristic r c (/f), superconductivity does not allow for SG order, while for smaller values a very 
detailed phase diagram arises with entangled spin glass and superconducting transitions. Our results 
for the Green's functions show that superconductivity occurring in the magnetic interaction band 
is of gapless type with a crossover from strongly gapless, within a certain range below T c , to very 
weakly gapless in a wide low temperature regime, and hardgapped at T — 0. 

PACS numbers: 71.27.+a, 74.25.Dw, 75.10.Nr 



Phase diagrams of HighTcSuperconductors (HTS) 
have been reported which show a spin glass phase in be- 
tween the antiferromagnetic and superconducting ones 
Strontium doped HTS are the most prominent 
examples, where antiferromagnetism gives way to clear 
signatures of spin glass (SG) like behaviour before su- 
perconductivity sets in at higher (hole)-doping. A typi- 
cal SG order parameter was identified at moderately low 
temperatures T < 8K Q and even an infiltration within 
the superconducting domain at lowest temperatures was 
described [0 . Classes of HTS exist as well, which do not 
seem to show spin glass and/or intermediate phase, or it 
might have been impossible to detect it up to now. Ex- 
istence and nature of intermediate spin glass or SG-alike 
phases in a certain doping range must be expected to be 
important features of strongly correlated systems. Ex- 
periments revealing close relations between magnetism 
and superconductivity in heavy fermion (HF) systems 
addressed coexistence, phase separation, and pairbreak- 
ing of local pairs by frozen moments for example . 
Most of recent theories for HTS materials focussed on the 
destruction of antiferromagnetism under doping. View- 
ing an intermediate phase from the superconducting side, 
the appealing concept of a nodal liquid [|7j was developed, 
leaving aside the notion and role of quenched disorder 
and randomness. Since we wish to deal here with su- 
perconductivity transitions under participation of a spin 
glass, a disorder model is a natural choice. Arguments on 
an important role of disorder in HTS and in HF-systems 
were provided experimentally and by theoretical reason- 
ing Aspects such as non Fermi liquid behaviour, 
seen to arise in the vicinity of spin glass order l| , also 
support this point of view. 

In this Letter we present detailed results for spin glass to 
superconductor transitions in a single-species fermionic 
model, which treats frustrated magnetic- and attractive 
interaction on the same footing. Unique features of the 



phase diagram are derived analytically and numerically. 
We shall observe that for certain interaction ratios the 
location of the spin glass bears resemblance to that of a 
logarithmic resistivity regime residing above a spin glass 
ordered phase at lower T in Sr-doped HTS Q. A fluc- 
tuational state of broken down spin glass order should 
contribute to transport properties seen in intermediate 
phases above Tf. In particular, SG-order was recently 
shown to affect transport properties strongly 0, an ef- 
fect that due to the random interaction can well have a 
weak localization precursor. 

We consider here a model described by the Hamilto- 
nian H = Hj v + Ht, where Hj v = ~^J2 Jij (J i< J j ~ 

E%4;4r a iT a .?l ~ P>Y, n h Tit = Ey CT V4 a i<x> and 
(7 = 7if— rii, a i°i and n = n T + n l denoting spin-, fermion- 
, and fermion-number operators respectively. The vari- 
ance J 2 of the frustrated and Gaussian-distributed mag- 
netic interaction Jij and its magnitude relative to that of 
the attractive coupling, v q =o / J, are relevant parameters 
of our analysis below, together with the chemical poten- 
tial and the related filling factor v(fi). We restrict the 
discussion to the small t/J regime, which implies that 
the selfconsistently determined magnetic band is much 
larger than the hopping bandwidth in general. Only deep 
within the superconducting regime, where the magnetic 
bandwidth almost shrank to zero, the single fermion hop- 
ping bandwidth becomes dominant. Local pairs can be 
delocalized due to finite range Vij and/or by arbitrarily 
weak fermion hopping ty . We employ a local pairing the- 
ory of superconductivity based on the order parameter 
A = (aii aii) and 011 a two-particle coherence length re- 
placing the usual one-particle length in the corresponding 
Ginzburg-Landau theory. One-fermion Green's functions 
JlOf are local as in the d — oo method Jll|,|l2| ■ Supercon- 
ductivity is of Bose condensation type and crosses over 
into BCS type behaviour (not considered here) in the 
limit of small ratios between v and hopping bandwidth. 
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The phase diagram illustrated in Figures 1 and 2 is ob- 
tained from the free energy for model TLj v given by 



f = \pj 2 {{q-l? 



-T 



00 dz 



1) 



2 / 2 lnC, 



(1) 



where C = cosh(/3ff) + cosh(/? v / |A| 2 + fi 2 ) exp(-±/?J 2 x) 
with x = 0(<1 — q) an d effective field H = J^/qz; the 
saddle point solutions q, q, and A are found from the 
coupled selfconsistent equations obtained by extremizing 
f. The physical solutions of these coupled equations are 
not easy to identify, even outside the SG-phase; multi- 
ple solutions exist and changeovers of stability occur as 
chemical potential /1, temperature T, or interaction ratio 
v/J are varied. This leads to the complexity of the phase 
diagram displayed in Figures 1 and 2. 
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FIG. 1. Phase diagram of competing spin glass (SG) and 
superconducting (SC)-transitions (thin lines) for various at- 
tractive couplings v (given by v — 2J/i e . p . at lines' right end- 
point (T = 0, fj, = fie.p.)), at fixed J = 1, and as a function of 
chemical potential fi (dark grey area: maximum SG-domain 
for v = 0). Thin lines enclose SC-phase which remove pieces 
of the (v — 0) SG-phase or suppress it totally at large enough 
v/J. The bold line I3 delimiting the light grey area joins tri- 
critical points (T c -maxima) on SC-critical curves locating 1st 
order transitions left of £3 

Figures 1 and 2 reveal unique features of the com- 
petition between SG- and SC-order. For example, an 
enhanced fermion concentration ^(/i) reduces the effec- 
tive spin density at larger and is seen to suppress the 
spin glass stronger than it eliminates the superconduct- 
ing phase as v(ijl) — ► or 2. Within a large region above 
the SG-phase, the SC-critical curves become deformed, 
as Fig.l shows for v/J < 4.5, due to the increasing 
spin glass fluctuations feeded by the random magnetic 
interaction. As the critical SC-curve passes through a 
maximum and starts to descend with decreasing /x, the 
SC-transitions change from second to first order. For 
still smaller ratios v/J the SC-line enters the v — SG- 
phase: in this case, magnetic moments freeze first and a 
discontinuous SG-SC transition follows at lower temper- 



ature (shown for v/J = 3.25,3.5,3.75). For still smaller 
v/J the 1st order SC-line falls rapidly to zero and the 
SG-phase prevents superconducting order up to a critical 
value fj, c (see v/J — 3,2.75,2.5). The replica symmetric 
solution also reveals the possibility of reentrant SC to 
SG-transitions, shown for v/J — 2.5, 2.75, 3 in Figure 1. 
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FIG. 2. Phase diagram close to the spin glass transition 
(bold line: pure SG-transition, bold-dashed line: paramag- 
netic stability limit) showing the stability limits of the super- 
conducting (SC) phase (dotted horizontal lines, given only for 
v/ J = 3.5, 3.0) and of the non-SC phase (dashed thin lines) bi- 
furcating at SC-tricritical points. SC-critical (continued thin) 
lines are shown down to small v/J. 



Figure 2 complements Figure 1 by adding, within a 
magnification of the SG-SC boundary, the stability lim- 
its of the superconductor and displays second order SC- 
transition lines crawling into the phase separation regime 
of the Ising spin glass for small v/J and T. For T — > 0, 
the quantum Parisi phase must be expected, which means 
that the paramagnetic stability limit should be shifted 
to the left. Lacking the full Parisi solution for /1 > we 
employ numerical solutions of fermionic TAP-equations 
Ip^fhil to arrive at a straight line estimate of the dis- 
continuous paramagnet-SG transition curve. Discontin- 
uous SC-SG transitions for j < 2.5 must occur in be- 
tween this curve and the SC-SG transition at j = 2.5. 
It is currently not possible to obtain a precise position of 
these discontinuous transitions at low T, since this would 
require to unite the dynamic mean field theory with 
the Parisi solution at infinite replica symmetry breaking 
(RSB) (despite the discontinuity, one step RSB is insuf- 
ficient at and near T = 0). An estimate can be obtained 
from the few percent rise of the SG-energy due to RSB 
]l5| . Furthermore, we arrive at the following conclusions: 

i) There is no coexistence of spin glass - and local pairing 
superconducting order parameter in zero magnetic field. 
The detailed analysis of the free energy and of all sta- 
bility conditions shows only SC-SG transitions between 
{q ^ 0, A = 0} and {q = 0, A ^ 0} for H = 0. 

ii) The transition between the two phases is always dis- 
continuous and exists only within a certain range of 
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chemical potentials (or filling). 

iii) For large enough v/J, like v/J > 4.13 at half-filling 
for example,the spin glass is prevented by the supercon- 
ducting transition, which is continuous for v/J> 4.55 
at half-filling; below this value it becomes discontinuous. 
The tricritical line, which separates these domains is in- 
cluded in Figure 1 for several values of v/J and as a 
function of the chemical potential. 

iv) Decreasing the temperature at fixed fi < .96 (see 
Figure 1), a second transition from SG to superconduc- 
tivity occurs at T c < Tf(v = 0). Figures 1 and 2 dis- 
play the exact numerical results of the replica symmet- 
ric theory. It is known that the spin glass free energy 
increases a few percent as an effect of replica symmetry 
breaking, which should slightly enlarge the superconduct- 
ing domain (these corrections increase with decreasing 
temperature). Recalling the change of the spin glass- 
ferromagnetic boundary under RSB [fl6[ the reentrant 
behaviour from superconductor to spin glass and back 
to superconductor may disappear as well, but the SC- 
SG boundary will not become a vertical line (unlike the 
boundaries between SG and ferro- or antiferromagnetic 
phase in a purely magnetic model [T^ ] ; 
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FIG. 3. Selfconsistent solutions for the superconducting or- 
der parameter A (continued lines), shown for v/ J = 3, q — 
and various /j, values. Unstable solutions (dashed lines, max- 
ima of f) explain stability limits and different discontinuous 
scenarios arising for sufficiently small fi. 



Figure 3 shows the effect of random magnetic fluctu- 
ations, described by q(T,fi), on the superconducting or- 
der parameter for exemplary values of \x and v/J = 3. 
The absence of coexistent order parameters allows to set 
q = 0. For [i = 0.8 a first order transition from A = 
to SC order occurs at a temperature below the tricriti- 
cal T3 = 0.475 J, the latter being independent of fi. For 
fx = J, the first order transition still occurs, but in addi- 
tion the non-SC phase becomes locally unstable between 
the two zeros of the finite A solutions. At low tempera- 
tures, A = is locally stable again, but the global min- 
imum remains at A 7^ 0. Beyond the tricritical point, 
(j, > fx 3 (j = 3) = 1.14 J in Figure 3, the stable SC solu- 
tion continuously decreases to A = 0. 



In a magnetic field new aspects arise: the transition tem- 
perature of the superconductor will be reduced and fi- 
nally vanish for H > H C 2, leaving a smeared spin glass 
transition for sufficiently small chemical potential. The 
overlap parameter q is nonzero in a field and then co- 
exists with A, since the field can penetrate the present 
type II superconductor partially. In this case the Almeida 
Thouless line can enter the superconductor, infiltrating 
crgodicity breaking there. 

We analyze the superconducting phase by means of the 
Green's functions: we derive the normal one G(e n ) , 
the anomalous one T{e n ), and relevant particle-hole- 
n p ^(w m ) and particle-particle bubble diagrams H pp (uj m ) 
for various limits. Superconductivity arising in the mag- 
netic band much larger than the hopping band is de- 
scribed by the (local) anomalous Green's function 

i A exp(/3 2 g/2)/(8gr 1 ) 



T{e) = 



cosh(/3ri) 
(e + iXr\ 



• exp(/? 2 q72) 



(2) 



{A,s} = ±l 

i A cosh(/3ri)/(4gri) v - 



1, 



3 (e + iXri + isf3q) 
2' 



^21 



2q 



cosh(/?n) 



■exp(^g/2)^ 



(e + iXn) U 



j 3 (j_ 

'2' 2q 



where n = ^//i 2 + A 2 , e = (2n + l)7rT, J = 1, U denotes 
the hypergeometric U-function (Kummer function) [p"7| , 
and the SC order parameter obeys A = vTJ2 e T(e). 





FIG. 4. Density of states p(E) (left) and TZe^iE)} (right), 
for (v — 5 J, fi — 0,T C ~ 1.1L7) and real energies E, show 
crossover from strongly to weakly gapless superconductivity 
for decreasing temperatures T/J = 1.1, 1.0, .9, .8 and .7 

The solution for G(e), being similarly given in terms of 
the Kummer function, yields the density of states plotted 
together with lZe{J-} in Figure |] (weak fermion hopping- 
and dynamic corrections from H p h(u)) and x(w) are neg- 
ligible here). These plots of course employ the stable 
selfconsistent solutions for q(T,ii) and A(T, //) inserted 
in (2). They demonstrate the crossover from strongly 
gapless superconductivity just below T c to a pronounced 
pseudogap. Although invisibly small for temperatures 
lower than roughly 20% below T c , the density of states 
remains nonzero in the pseudogap regime and vanishes 
there only at T — 0. The term gapless superconduc- 
tor should be used with care, since over a wide temper- 
ature range the gap looks almost perfect. The round- 
ing of the superconducting gap edges and the retarded 
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way it opens up below T c is reminiscent of the magnetic 
gap found below spin glass transitions ||. As the super- 
conducting gap sharpens at lower T, the magnetic band 
narrows and eventually the fermion hopping (bandwidth) 
starts to dominate. Then the anomalous Green's func- 
tion crosses over into the hopping band solution which 
we derived for gaussian random tij as 



A 



2w 2 y/2w 2 



A yjw 2 + e 2 



r(e) 



Ve 2 + A 2 



(3) 



where r(e) = \J Aw 2 it 



with w 2 = 4(t?-). Equation 



+ A 2 ) + {-w 2 + e 2 + A 2 + /i 2 ) 2 
, , .. ^.j .^vion (^) is an exact local solution 
which holds for arbitrary bandwidth w and chemical po- 
tential (i, extending an earlier result for half-filling [ fujl . 
The transition temperature derived from (^) illustrates 
in Figure | the crossover from linear Bose condensation 
T c (w)-dependence to exponential BCS-type. 




FIG. 5. Complete crossover from Bose condensation type 
to BCS behaviour illustrated in 3D plot of T c (v) as w/v in- 
creases; typical values for both linear and exponential depen- 
dence are given as a function of A = vpF 

The Bose-BCS crossover can be followed through the 
whole parameter range due to the local property of the 
one-particle Green's functions, a feature appreciated as 
well in the d = oo method for clean systems |ll],[l^] • The 
suppression of one particle phase coherence, due to infi- 
nite dimensions in clean systems or by symmetry require- 
ment in the quenched average used here, results in type 
II superconductivity with two particle coherence length 
replacing the standard one in the Ginzburg Landau the- 
ory (many details including paramagnetic pairbreaking 
were published in The puzzling coexistence phe- 

nomena in zero field, absence of (q ^ 0, A ^ 0) so- 
lutions but presence of phase separation regimes, and 
the related discontinuous transition between the mag- 
netic and superconducting phases is further elucidated 
by the bubble diagram H PP (T, lo). At half- filling for 
example, the condition 1 = vH pp (T c ,lu = 0) yields 
T c = (1 - q(T c ))l , (q = T X | w=0 ). Equating this with 
the spin glass freezing temperature Tf = Jq(Tf), ie solv- 
ing formally T c — Tf to search for simultaneous onset 
of both types of order, results in q = v/(v + J). Since 
q(Tf) = .6767 one finds v = 4.1862 (let J = 1). Since 
max{n(T,^ = 0,uj = 0)} = .242 at T = .731 no solu- 
tion exists for v < i> m i n = 4.14, while two solutions are 
found for v = 4.186 > i> m i n , only those with T > .731 be- 
ing free energy minima. At half-filling, the thermal first 



order superconducting transition occurs at temperature 
T c i > Tf for v > 3.8, but the stability limit of the normal 
phase is absent for 3.8 < v < 4.14; spin glass order will 
not reappear down to lowest temperatures. For v < 3.8 
the spin glass transition occurs at Tf = .6767 with a first 
order SG-SC transition to follow at lower temperature. 
For the sake of transparency we discussed only fully frus- 
trated magnetic interactions and zero field phenomena. 
Subsequent antiferromagnetic-spin glass-superconductor 
transitions, as \x increases from zero are not simply ob- 
tained once the Jij is given an antiferromagnetic part. 
Thus, models extended by the Hubbard interaction, but 
also the allowance for different symmetries of the order 
parameter, and the dynamic quantum Parisi phase are 
examples for future research on links between antiferro- 
magnetism, spin glass order, and superconductivity. 
The model we analyzed here proved to have a special 
type of gapless superconductivity due to the vicinity of 
spin glass order and due to correlations induced by the 
frustrated magnetic interaction. 
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[1] 

[2] 
[3] 
[4] 
[5] 

[6] 



F. C. Chou, N. R. Belk, M. A. Kastner, R. J. Birgeneau, 

and A. Aharony, Phys. Rev. Lett. 75, 2204 (1995) 

D. J. Scalapino, Phys. Rep. 250, 329 (1995) 

W. Brenig, Phys. Rep. 251, 155 (1995) 

J. A. Mydosh, J. Magn. Magn. Mat. 157/158, 606 (1996) 

I. Korenblit, V. Cherepanov, A. Aharony, and O. Entin 

- Wohlmann, |sond-mat/9709056| 

H. Spille, M. Winkelmann, U. Ahlheim, C. Bredl, F. 
Steglich, P. Haen, J. Migot, J. Thoulence, and R. 
Tournier, J. Magn. Magn. Mat. 76/77, 539 (1988)and 
references therein 



[7] 
[S] 

[9] 
[10] 

[11] 

[12] 

[13] 

[14] 
[15] 
[16] 
[17] 

[18] 



L. Balents, M.P.A. Fisher, C. Nayak, cond-mat/9803089 



A. Georges and A. Sengupta, Phys. Rev. B52, 10295 
(1995) 

R.Oppermann and B. Rosenow, Phys. Rev. Lett. 80, 4767 

(1998) and to appear in Phys. Rev. B (1998) 

R. Oppermann, Physica A167, 301 (1990) and references 

therein 

D. Vollhardt in Correlated electron systems, World Sci- 
entific, 57 (1994) 

A. Georges, G. Kotliar, W. Krauth, and M. Rozenberg, 
Rev. Mod. Phys. 68 13 (1996). 

D.J. Thouless, P.W. Anderson, and R.G. Palmer, Phil. 
Mag.35, 593 (1977) 



M. Rehker and R. Oppermann, :ond-mat/9806092 
G. Parisi, J. Phys. A13, 1101 (1980) 
K. Binder, A.P. Young, Rev. Mod. Phys. 58, 801 (1986) 
M. Abramovitz and I. Stegun, Handbook of Mathemati- 
cal Functions, Dover Publ. Inc, NY (1965) 
B. Rosenow and R. Oppermann, Phys. Rev. Lett. 77 
1608 (1996) 



4 



